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<D : Abstract 

Magnetization dynamics in the system of magnetic nanoclusters with randomly 
, distributed properties are studied by means of computer simulations. The main at- 

tention is paid to the possibility of coherent magnetization reversal from a strongly 
' nonequilibrium state with a mean cluster magnetization directed opposite to an ex- 

■ ternal magnetic field. Magnetic nanoclusters are known to be characterized by large 

2 I magnetic anisotropy and strong dipole interactions. It is also impossible to produce a 

number of nanoclusters with identical properties. As a result, any realistic system of 
nanoclusters is composed of the clusters with randomly varying anisotropics, effective 
^ ', spins, and dipole interactions. Despite this randomness, it is possible to find condi- 

CTn I tions when the cluster spins move coherently and display fast magnetization reversal 

due to the feedback action of resonator. By analyzing the influence of different cluster 



l/^ , parameters, we find their optimal values providing fast magnetization reversal. 

o 
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1 Introduction 



Magnetic nanoclusters are the objects enjoying rich and interesting properties with a variety 
of apphcations, as can be inferred form the review articles |TH5]. For example, they are 
used in magnetic chemistry as catalysts; in biomedical imaging for magnetic reading by 
magnetometers that measure the magnetic field change induced by nanoclusters; in medical 
treatment by employing alternating magnetic fields forcing the oscillation of the nanocluster 
magnetic moment that produces heat destroying ill cells or bacteria; in genetic engineering 
by attaching a nanocluster to a particular piece of the molecule and then removing it together 
with this piece; in waste cleaning by stacking nanoclusters to waste and then removing them 
together with waste; in information storage, processing, and quantum computing; and so 
on [iHS]. This is why the methods of governing magnetization dynamics of nanoclusters are 
so much important. 

The use of nano sizes of clusters is principal, since such clusters behave as large mon- 
odomain particles with a given large total spin. Typical radii of nanoclusters are between 1 
and 100 nm, containing from 100 to 10^ atoms. Respectively, the effective total magnetiza- 
tion of a cluster can be of order 100 to 10^ hb- The maximal size of a cluster, when it is in a 
monodomain state, such that the spins of atoms, composing the cluster, sum up coherently, 
forming an effective total spin, is characterized by the coherence radius. The clusters, whose 
sizes are larger than the coherence radius, decompose into several domains with oppositely 
directed magnetization, so that the total cluster spin is zero. 

Employing nanoclusters for the purpose of information storage, processing, and quantum 
computing meets the requirements that contradict each other. From one side, for information 
storage, it is required that the cluster magnetization could be well frozen in a given direction, 
which needs the existence of a sufficiently large magnetic anisotropy. But, from another 
side, for information processing, it is necessary to be able to quickly reverse the cluster 
magnetization, which is hindered by this anisotropy. 

Nanoclusters enjoy the property to get their magnetization frozen at temperatures lower 
than the blocking temperature that is of order of 10 - 100 K. Then to reverse the cluster 
magnetization, one has to overcome the magnetic anisotropy. So, for information storage, 
one needs large magnetic anisotropy, while for information processing, the anisotropy has 
to be suppressed to achieve fast magnetization reversal. To accomplish such a reversal, 
one applies external transverse fields pushing the cluster magnetization. The discussion of 
different methods of magnetization reversal by means of external magnetic fields can be 
found in Refs. [2H6|. 

A very efficient method for realizing fast magnetization reversal of nanoclusters is based 
on the use of feedback field from a resonant electric circuit coupled to the ensemble of nan- 
oclusters [BHTU]. Actually, the idea that a resonant electric circuit can drastically shorten 
the relaxation time of a spin system was advanced by Purcell [llj, who illustrated it for an 
ensemble of nuclear spins. This mechanism has later been considered by Blombergen and 
Pound [12]. The influence of the coupling of a magnetic sample with a resonant circuit on 
magnetization dynamics is called the Purcell effect. One sometimes calls it the radiation 
damping, following Bloembergen and Pound [12j. However the latter term is rather a mis- 
nomer, as has been stressed by a number of researchers [TBHIS], since the coupling of a spin 
system with a resonator does strongly influences spin dynamics, but not necessarily damping 
it. For instance, this effect enhances nuclear magnetic resonance signals, which is widely used 
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in NMR techniques |16I - [T8] . The term "radiation damping" is also confusing because of mak- 
ing impression that these are spins themselves that cooperate by interacting with each other 
through a common radiation field. Such a radiation correlation between radiating atoms is 
the essence of the Dicke effect [19]. But for magnetic particles, as has already been stressed 
by Purcell [11], such a radiation collectivization of spin motion is absolutely negligible. In 
atomic physics, one clearly distinguishes these two principally different physical effects. And 
the term Purcell effect is employed for describing the influence of a cavity resonator on atomic 
radiation, which is the main part of the cavity quantum electrodynamics [20] - [22j . 

The role of the Purcell effect on the spin dynamics of nuclear spins is well studied. In 
resonance experiments, it enhances the NMR signals [T614T8] . For strongly nonequilibrium 
systems of polarized nuclei, it leads to fast magnetization reversal [23H27] that has been dis- 
covered in experiments [23] and later confirmed in other experimental studies (see references 
in the review article [1]). 

Magnetization dynamics in an ensemble of nanoclusters is essentially different from that 
of nuclear systems. The basic differences are as follows. 

(i) Nanoclusters possess rather large magnetic anisotropy that hinders the possibility of 
simple regulation of spin motion. 

(ii) Having large effective spins, nanoclusters also have strong spin dipole interactions, 
which results in short dephasing time. 

(iii) The most important difficulty for organizing collective spin motion in a system of 
nanoclusters is the problem of cluster inhomogeneity, since nanoclusters, being prepared by 
any of the known methods, whether by thermal decomposition, or microemulsion reactions, 
or by thermal spraying, are not identical particles, as nuclei would be. But nanoclusters 
differ by their shapes and sizes, which results in the difference in the value of their spins, 
dipole interactions, and of their anisotropics. 

It is the aim of the present paper to investigate the magnetization dynamics in a realistic 
inhomogeneous ensemble of nanoclusters exhibiting all their typical properties of anisotropy 
and dipole interactions. Since analytic investigation of such an inhomogeneous system is too 
much complicated, we resort to computer modelling. Because the problem of magnetization 
reversal is of special interest for many applications, such as information recording and pro- 
cessing, we pay the main attention to finding the optimal conditions, when the magnetization 
reversal is fast and as close to complete as possible. 

2 Realistic model of nanocluster system 

The sample formed by a system of nanoclusters is described by the Hamiltonian 



where the first term characterizes single nanoclusters and the second term describes their 
dipole interactions. 

The typical Hamiltonian of a nanocluster has the form 
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in which the first term is the Zeeman energy, and other terms describe the energy due to 
magnetic anisotropy, with the corresponding anisotropy parameters [IHS]- 
The interaction Hamiltonian is caused by dipole spin interactions 



Hij = ^ DifS!;'Sj , (3) 

a/3 



with the dipolar tensor 



{Sap - 3<,.ng) , (4) 



where 



II ^ij 

^ij — I ^ij I ; ^ij — ; ^ij — ■ 



The total magnetic field, acting on the system, is the sum 

B = Boe, + He, (5) 

of an external constant field Bq and the resonator feedback field H. 

The considered sample is inserted into a coil of an electric circuit with an attenuation 7 
and natural frequency u. Moving spins of the sample create electric current in the coil that, 
in turn, produces the feedback magnetic field acting on these spins. The equation for the 
feedback field follows from the Kirchhoff equation and can be written jH[26l[27] as 



^ + 2lH + co' 1^ H{t') dt' = -Anv ^ , (6) 

where rj is filling factor and the electromotive force is caused by the moving average magne- 
tization density 

j 

The equations of spin dynamics are obtained in the following traditional way [Hll4y 28 [ l29]. 
First, we write the Heisenberg equations of motion for the spin components S"", in which 
the index j = 1,2, . . . , N enumerates nanoclusters and a = x,y, z. Then these equations are 
averaged using semiclassical approximation, and spin attenuation is taken into account by 
means of the second-order perturbation theory, as is described in full details by Abragam and 
Goldman [lll[29]. This procedure results in the equations that can effectively be obtained 
from the evolution equations 

zhj^S^ = [S^,H]-zT2S^, 



th-S] = [S],H]~tT2S] , 

thj^S] = [S^,H]-t^,{S]-S) , (8) 

with the straightforward use of the semiclassical approximation |ll[Tll|28l[29]. Here r2 is 
a transverse attenuation parameter, 71 is longitudinal attenuation parameter, and S = 
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Ylj'=i^j/-^ is average spin value. In the case of strong initial polarization, the transverse 
attenuation is characterized [ll[Tlll271[2n] by the expression 

r2 = 72(l-s'), (9) 



where s is the reduced average spin 



TV 



(10) 



j=i ■ J 



and the natural width is 



72 - ^ = Phils . (11) 

In the following numerical calculations, we shall use the reduced dimensionless attenua- 
tion parameters, measured in units of the Zeeman frequency 

^„ ^ _ = \^^\Bo , ifio = -2fiB = -hhs\) . (12) 

n 

Our goal is to investigate spin dynamics of nanoclusters with realistic parameters. Thus, 
the values of the anisotropy parameters, typical for nanoclusters, such as Co, Fe, and Ni 
nanoclusters, are 

-^ = lo-^ #^ = lo-^ ^ = 10-^0. (13) 

A172 ri'j2 '172 

Studying spin dynamics, we shall take into account that real nanoclusters are not com- 
pletely identical with each other, but exhibit the dispersion in their anisotropy parameters 
and spin values. The main aim of our study is to analyze how the parameter dispersion 
influences spin dynamics and to find conditions allowing for effective spin reversal. 



3 Introduction of dimensionless system parameters 

For numerical analysis, it is convenient to pass to dimensionless quantities. To this end, we 
define the dimensionless resonator feedback field 

VH^^ (14) 

and the reduced transverse attenuation, caused by dipole interactions, 

P. = - . (15) 

Wo 

Also, we introduce the dimensionless anisotropy parameters 

Pa = — , Pb = — , Pc^ — , (16) 

Uq 

expressed through the anisotropy frequencies 

0JA = 2S^, ujB = 2S^, ujc = 2S'^. (17) 
n n n 
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In the equations of motion, we employ the reduced spin, with the components 

j=l 

Since the nanocluster spins Sj are large, of order 10^, the semiclassical approximation is well 
justified. This allows us to treat e as a classical vector. The evolution will be considered 
with respect to the dimensionless time 

t = uot . (19) 

At the initial moment of time, there is no feedback field, which assumes the initial 
conditions 

Ph{0) = , Ph{0) = , (20) 

where the overdot implies time derivative. 

The initial conditions for the spin polarization e^(0) are constructed as follows [301 l3T] . 
For the given value of an initial polarization 6^(0), a variety of the initial orientations for 
the individual vectors Sj{0) are admissible. These orientations can be prescribed by a kind 
of Monte Carlo techniques. First, a random configuration of vectors Sj is taken and the 
corresponding total polarization is evaluated. A new direction is chosen randomly for each 
spin, and the new total polarization is calculated. If it is less than the initial one, the 
array with the changed magnetic moment direction distribution is chosen as the second 
iteration, otherwise, it is rejected. This procedure is repeated until the system achieves the 
required average polarization ez{0) playing the role of the initial condition. The initial spin 
polarization is assumed to be directed opposite to the external magnetic field Bq. 



4 Optimal conditions for spin reversal 

We accomplish numerical solution of the evolution equations for = 15^ = 3375 spins, 
searching for the conditions of the optimal spin reversal. The Zeeman frequency ojq is taken 
to be in resonance with the circuit natural frequency u. 

The peculiarity of spin reversal depends on the system parameters. For instance, the op- 
timal value of the resonator attenuation 7, defined by the resonator circuit quality, essentially 
depends on the given transverse attenuation 72, caused by spin dipole interactions. 

Let us fix 72 = O.OlcJo, which yields pd = 0.01, and let us take the anisotropy parameters 
typical for Co, Fe, and Ni nanoclusters, as defined above. Below the blocking temperature, 
the longitudinal relaxation is suppressed, so that 71 is much smaller than 72, which allows 
us to set 71/72 = 10~^. The coil filling factor is assumed to be close to one. The nanocluster 
spins are randomly distributed by the normal law with mean 5* = 1200 and variance 6S = 
400. Dynamics of spin polarization, e^, is shown in Fig. 1 for different resonator attenuations 
7. Varying 7, we are looking for its optimal value that provides the maximal permanent 
reversal, without superimposed oscillations. It is this regime that is optimal for fast and 
stable information processing [22] • As is seen, for the given ratio 72/wo; the optimal value 
of the resonator attenuation is 7 ~ 0.35ci;o. Below this value, the reversal is slower and the 
reversed spin value is smaller. For very low 7, there is no reversal at all. Above the optimal 
value of 7, there arise oscillations, and the reversed spin value diminishes. 
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The spin reversal is caused by the resonator feedback field. This is ilhistratcd in Fig. 2, 
where it is seen that the reversal occurs when the amplitude of the fccdback-ficld oscillations 
is maximal and almost reaches the strength of the external field Bq. For each given ratio 
Pd = li/'^O) the corresponding optimal value of 7 is chosen. The stronger external field Bq, 
that is, the larger uiq, hence, smaller p^, the more effective is the reversal, though the reversal 
time increases. Spin magnitudes are randomly distributed, as in Fig. 1. 

The dependence of the optimal 7* and the related value of the reversed spin s* = e*, as 
functions of the ratio pd = 12/^0, are presented in Fig. 3. The larger the ratio pd, the larger 
the optimal 7*, but the smaller the value of the reversed spin s*. 

In addition to spin dispersion, there can exist the dispersion of the anisotropy parameters. 
We study spin dynamics, when the spins as well as all anisotropy parameters are distributed 
by the normal law. The spin mean is taken as S* = 1200 and variance as 6S = 960. The 
mean values of the anisotropy parameters are assumed to be typical for the Co, Fe, and Ni 
nanoclusters, as in Eq. (13), with the relative variance 

SD _ 5D2 _ _ Q ^ 
D D2 D4 

We accomplish 30 realizations of spin dynamics varying the anisotropy parameters, with ran- 
dom spin distribution in each of the realizations. The results, averaged over the realizations, 
are shown in Fig. 4. We see that a rather strong dispersion of the system parameters does 
not preclude spin reversal. 

To understand when the spin reversal could be blocked by anisotropy, we vary the pa- 
rameters of the latter in a wide range, looking for the values at which the reversal becomes 
blocked. Varying one of the anisotropy parameters, we keep others fixed. Spins are randomly 
distributed with mean S = 1200 and relative variance 0.3. Fig. 5 illustrates the results, from 
which the blocking anisotropy parameters are evaluated as 

^ ^ 10-2 , ^ ^ 0.5 X 10-^ , ^ ~ 10-« . 
ft72 n.72 ^72 

These blocking parameters are essentially larger than the typical values in Eq. (13). Hence, 
the typical anisotropy does not block spin reversal. 

With increasing external field Bq, that is, diminishing the ratio pd = 72/(^0, the blocking 
anisotropy parameters decrease. This is illustrated by Fig. 6 that is to be compared with 
Fig. 5. So, increasing the external field suppresses the role of the anisotropy. 

Finally, we study whether a strong spin dispersion can destroy the coherent spin dynam- 
ics and infiuence spin reversal. For this purpose, we consider a random distribution of spins, 
with mean S = 1200 and very large relative variance 6S/S = 1. The results of numerical 
simulations are given in Fig. 7 for the ratio pd — 72/1^0 = 0.01, with the corresponding 
optimal j/ujq — 0.3, and varying anisotropy parameters. The results demonstrate that, for 
small anisotropy parameters, spin dispersion plays practically no role. This role increases for 
larger values of the anisotropy parameters. This fact finds straightforward explanation from 
the structure of the spin equations of motion. When the anisotropy parameters are much 
smaller than the Zeeman frequency, the spin rotation is governed by the same cuq depending 
only on the external field Bq, but independent of the spin lengths. But when the anisotropy 
parameters are sufficiently large, approaching their blocking values, then the anisotropy 
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disturbs the effective rotation frequencies, so that different spins rotate with different fre- 
quencies. This is equivalent to the appearance of inhomogeneous broadening. Fortunately, 
the typical nanoclustcr anisotropy parameters are much lower than their blocking values. 
So, for a sufficiently strong external field, all spins rotate with almost the same frequency 
ujq and their dispersion even being quite large, does not much influence their motion. Then 
the coherent spin dynamics can be realized, resulting in a fast spin reversal. 

5 Conclusion 

We have studied the magnetization dynamics in an ensemble of magnetic nanoclusters, keep- 
ing in mind a realistic situation, when the nanoclusters are characterized by a microscopic 
Hamiltonian, with well defined parameters. Another peculiarity of realistic nanocluster en- 
sembles is the dispersion in the values of their spins and anisotropy parameters, which we 
also take into account. We have accomplished a series of computer simulations varying the 
system properties and searching for the conditions when the magnetization reversal is opti- 
mal, in the sense of being fast, maximal, and quasi-stationary, without oscillations. Such a 
process of magnetization reversal is necessary in a variety of applications. For instance, for 
realizing effective quantum information processing. 

For numerical values, wc have chosen the nanoclustcr parameters typical for Co, Fe, and 
Ni nanoclusters. It turns out that if the external magnetic field defines the Zeeman frequency 
that is essentially larger than the effective anisotropy frequencies, then the dispersion of spin 
magnitudes plays practically no role. This also concerns the dispersion in the values of the 
anisotropy parameters. The external magnetic field, sufficient for this effect is 1 T, which 
corresponds to the Zeeman frequency uq ~ lO^-*^ Hz. 

We have shown that the magnetization reversal is caused by the resonator feedback field 
that is self-organized by moving spins. The value of the feedback field, in the moment of 
spin reversal, almost reaches the value of the external field. Fast magnetization reversal 
can be achieved solely by moving spins themselves, without involving strong transverse 
magnetic fields. For such nanoclusters as Co, Fe, and Ni, the reversal time can be of order 
10^^^ s. This conclusion suggests a convenient mechanism for the efficient manipulation of 
nanoclustcr magnetization, which is of high importance for a variety of applications. 
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Figure Captions 



Fig. 1. Dynamics of tlie average spin polarization Cz = s, for the ratio 72/co'o = 0.01, 
as a function of dimensionless time, for random spins with the mean S — 1200 and relative 
spin variance 5S/S = 0.333 for different values of the resonator attenuation: (a) 7/wo = 
(dashed-dotted line), 7/^0 = 0.04 (dashed line), 'j/uq = 0.12 (solid line); (b) ^/ojq — 0.3 
(dashed- dotted line), ■y/cuo = 0.4 (dashed line), 7/a;o = 0.8 (solid line). 

Fig. 2. Average spin polarization = s (dashed line) and feedback field pn (solid line) 
as functions of dimensionless time, for different ratios = 72/1^0, with the corresponding 
optimal 7: (a) pd = 0.1,7/wo = 0.7; (b) pd = 0.01,7/a;o = 0.3; (c) Pd = 0.001, 7/a'o = 0.1. 
Spins are randomly distributed as in Fig.l. 

Fig. 3. Dependence of the optimal 7* and the related value of the reversed spin s* — e*, 
as functions of the ratio pd = 72/(^0 • 

Fig. 4. Influence of the combined dispersion of spins and anisotropy parameters on the 
dynamics of spin reversal for pd = 72/1^0 = 0.01, with the corresponding optimal 7/^0 = 0.3: 
(a) spin polarization s = Cz (dashed line) and feedback field pn (solid line) as functions 
of the dimensionless time for random spins with the relative spin variance 6S/S — 0.8 and 
fixed anisotropy parameters; (b) Spin polarization s — Cz, averaged over 30 realizations of the 
anisotropy distribution with the relative variance 0.1, as a function of the dimensionless time; 
(c) feedback field as a function of the dimensionless time, averaged over 30 realizations 
of the anisotropy distribution. 

Fig. 5. Spin polarization s = e^;, as a function of dimensionless time, for the ratio 
Pd = 12/^0 = 0.01, with the corresponding optimal 7/1^0 = 0.3, for varying anisotropy 
parameters: (a) pa = 0.1 (sohd line), pA = 0.2 (dashed line), pa = 0.3 (dashed-dotted line). 
Pa — 0.4 (dotted line); (b) ps = 0.5 (solid line), ps = I (dashed line), ps = 2 (dashed- 
dotted line), Pb — ^ (dotted hne); (c) pc — 0.1 (sohd hne), pc — 0.2 (dashed line), pc — 0.3 
(dashed-dotted line), pc = 0.4 (dotted hne). In each case, spins are randomly distributed 
with mean 1200 and relative variance 0.3. 



Fig. 6. Spin polarization s = e^, as a function of dimensionless time, for the ratio 
Pd — 72/^0 = 0.01, with the corresponding optimal j/ujq = 0.3, for varying anisotropy 
parameters: (a) pa = 0.01 (solid line), pa = 0.05 (dashed line), pa = 0.1 (dashed-dotted 
line). Pa = 0.2 (dotted line); (b) ps = 0.05 (sohd line), = 0.1 (dashed line), ps = 0.2 
(dashed-dotted line), pb = 0.5 (dotted hne); (c) pc — 0.01 (sohd line), pc — 0.05 (dashed 
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line), PC = 0.1 (dashed-dotted line), pc = 0.2 (dotted line). In each case, spins are randomly 
distributed with mean 1200 and relative variance 0.3. 



Fig. 7. Spin polarization s = as a function of dimensionless time, for the ratio 

Pd = 12/^0 = 0.01, with the related optimal j/ujq = 0.3, for varying anisotropy parameters. 
Solid line describes the case of strong spin dispersion, with mean spin S = 1200 and relative 
variance 6S/S = 1. Dashed line corresponds to the case without spin dispersion. The 
anisotropy parameters are: (a) pa = 0.01, pb = 0.01, pc = 0.001, (b) pa = 0.05, ps = 0.1, 
PC = 0.05, (c) PA = 0.1, Pb = 0.2, pc = 0.1, (d) pa = 0.2, pb = 0.7, pc = 0.2, (e) pa = 0.05, 
PB = 0.01, PC = 0.001, (f) PA = 0.01, PB = 0.1, PC = 0.001, (g) PA = 0.01, PB = 0.01, 
PC = 0.02. 
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Figure 1: Dynamics of the average spin polarization = s, for the ratio 72/^0 = 0.01, as 
a function of dimensionless time, for random spins with the mean 5* = 1200 and relative 
spin variance 5S/S = 0.333 for different values of the resonator attenuation: (a) 7/070 = 
(dashed-dotted line), ^/uq — 0.04 (dashed line), ^/loq — 0.12 (solid line); (b) ^/ujq — 0.3 
(dashed-dotted line), ^/uo — 0.4 (dashed line), ^/loq — 0.8 (sohd line). 
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Figure 2: Average spin polarization = s (dashed line) and feedback field p^j (solid line) 
as functions of dimensionless time, for different ratios pd = li/^O) with the corresponding 
optimal 7: (a) Pd = 0.1,7/(^0 = 0.7; (b) pa = 0.01, 7/(^0 = 0.3; (c) pa = 0.001, 7/^0 = 0.1. 
Spins are randomly distributed as in Fig.l. 
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Figure 3: Dependence of the optimal 7* and the related value of the reversed spin s* = e*, 
as functions of the ratio pd = 72/i^o- 
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Figure 4: Influence of the combined dispersion of spins and anisotropy parameters on the 
dynamics of spin reversal for = 12/^0 = 0.01, with the corresponding optimal 7/(^0 = 0.3: 
(a) spin polarization s = (dashed line) and feedback field pn (solid line) as functions 
of the dimensionless time for random spins with the relative spin variance 5S/S — 0.8 and 
fixed anisotropy parameters; (b) Spin polarization s = e^, averaged over 30 realizations of the 
anisotropy distribution with the relative variance 0.1, as a function of the dimensionless time; 
(c) feedback field pn as a function of the dimensionless time, averaged over 30 realizations 
of the anisotropy distribution. 



16 




Figure 5: Spin polarization s = e^, as a function of dimensionless time, for the ratio pd = 
72/'^o = 0.01, with the corresponding optimal ■y/coo = 0.3, for varying anisotropy parameters: 
(a) Pa = 0.1 (solid line), pa = 0.2 (dashed line), pa = 0.3 (dashed-dotted line), pa = 0.4 
(dotted hne); (b) pb — 0.5 (sohd hne), pb — 1 (dashed line), pb — 2 (dashed-dotted line), 
Pb = 3 (dotted line); (c) pc = 0.1 (solid line), pc — 0.2 (dashed line), pc = 0.3 (dashed- 
dotted line), Pc = 0.4 (dotted line). In each case, spins are randomly distributed with mean 
1200 and relative variance 0.3. 
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Figure 6: Spin polarization s = e^, as a function of dimensionless time, for the ratio pd = 
72/'^o = 0.01, with the corresponding optimal ""y/oJo = 0.3, for varying anisotropy parameters: 
(a) Pa = 0.01 (solid line), pa = 0.05 (dashed line), = 0.1 (dashed-dotted line), pa = 0.2 
(dotted line); (b) pb — 0.05 (sohd line), pb — 0.1 (dashed hne), pb — 0.2 (dashed-dotted 
line), Pb = 0.5 (dotted line); (c) pc = 0.01 (solid line), pc = 0.05 (dashed line), pc = 0.1 
(dashed-dotted line), pc = 0.2 (dotted hne). In each case, spins are randomly distributed 
with mean 1200 and relative variance 0.3. 
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Figure 7: Spin polarization s = e^, as a function of dimensionless time, for the ratio = 
72/(^0 = 0.01, with the related optimal 7/a;o = 0.3, for varying anisotropy parameters. Solid 
line describes the case of strong spin dispersion, with mean spin 5" = 1200 and relative 
variance bSjS — 1. Dashed line corresponds to the case without spin dispersion. The 
anisotropy parameters are: (a) 'Pa = 0.01, pb = 0.01, pc = 0.001, (b) pa = 0.05, pb = 0.1, 
PC = 0.05, (c) PA = 0.1, Pb = 0.2, pc = 0.1, (d) pa = 0.2, pb = 0.7, pc = 0.2, (c) pa = 0.05, 
PB = 0.01, PC = 0.001, (f) PA = 0.01, PB = 0.1, PC = 0.001, (g) PA = 0.01, PB = 0.01, 
PC = 0.02. 
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